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ABSTRACT 



We present a multiwavelength study of the black hole X-ray binary V404 Cyg 
in quiescence, focusing upon the spectral energy distribution (SED). Radio, optical, 
UV, and X-ray coverage is simultaneous. We supplement the SED with additional 
non-simultaneous data in the optical through infrared where necessary. The compiled 
SED is the most complete available for this, the X-ray and radio brightest quiescent 
black hole system. We find no need for a substantial contribution from accretion light 
from the near-UV to the near-IR, and in particular the weak UV emission constrains 
published spectral models for V404 Cyg. We confirm that no plausible companion 
spectrum and interstellar extinction can fully explain the mid-IR, however, and an IR 
excess from a jet or cool disc appears to be required. The X-ray spectrum is consistent 
with aT~2 power-law as found by all other studies to date. There is no evidence for 
any variation in the hardness over a range of a factor of 10 in luminosity. The radio 
flux is consistent with a flat spectrum (in /„). The break frequency between a flat and 
optically thin spectrum most likely occurs in the mid or far-IR, but is not strongly 
constrained by these data. We find the radio to be substantially variable but with no 
clear correlation with X-ray variability. 
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1 1 INTRODUCTION 

The question of how accretion onto a black hole occurs has 
been of great and ongoing interest. Observations of accreting 
sources provide a test of astrophysics in the strong gravity 
regime and will potentially yield a signature of the presence 
of an event horizon around black holes. Recent advances 
in X-ray technology have revealed that accretion can oc- 
cur at extremely low rates across the range of black hole 
masses from stellar remnants to supermassive black holes 
at the centre of our own and other galaxies. However, the 
nature of that low-level accre tion remains unclear (see e.g., 
iNaravan fc McClfntockll200^) . Accretion flows at low lumi- 
nosities do appear to be radiatively inefficient, but we have 
yet to determine to what extent that arises from advec- 
tion of hot gas through the event hor izon, and what rol e 
jets may play in carrying away energy l|Fender et ahl feoOS). 
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These low-level accreting black holes emit across the electro- 
magnetic spectrum from radio to X-rays, so multiwavelength 
studies can be used to disentangle different sources from dif- 
ferent regions of the inflow, outflow, or jet, and to establish 
causal connections between behavior in these regions. 

Among the stellar mass black hole population, the 
best studied quiescent system for this purpose has been 
V404 Cyg, which contains a compact object with a mass 
function of 6.08 ± 0.06 M Q (ICasares fc Charles! ll994T ) - one 
of the most secure stellar mass black hole mass determina- 
tions we have. V404 Cyg has the longest orbital period of 
the black hole X-ray tran sients (BHXRTs), at 6.5 days, with 
a K0IV companion star l|Casares et al.|[l993l ). It is also the 
most luminous stellar mass black hole in quiescence with 
L x ~ 10 33 - 10 34 ergs -1 , ~ 1000 x brighter tha n the proto- 
typical BHXRT, A 0620-00 l|Garcia et all 120011 ). These ex- 
tremes are unlikely to be coincidental and the high Lx may 
be a consequence of the long orbital period, the evolved 
donor star, and/or the high mass of the black hole. 

Attempts to relate observations of quiescent black hole 
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binaries to theoretical models of accretion flows have fo- 
cused upon fitting the broad-band spectral energy distri- 
bution, ideally with coverage from radio to X-ray energies. 
V404 Cyg has been a prime target due to its X-ray bright - 
ness (|Naravan et al.lll996l . [l997l ; lQuataert fc Naravanll 19991) . 
These studies mainly focused on pure accretion models and 
only considered the optical to X-ray regime, as radio data 
were limited or non-existent. Attempts to appl y advective 
mode ls to the more luminous hard state (e.g. Esin et ail 
2001) brought home the fact that these models could not 
explain the ubiquitous radio emission in this state, that was 
sometimes seen to extend as far as the near-IR band. A more 
complete description of the low-state SED a ppears to re- 
quire a coupled accretio n plus outflow model (iMarkoff et all 
l200ll ; lYuan et al.|[2005l ). More recent studies have revealed 
significant radio e mission in quiescence to o, from V404 Cyg 
(|Gallo et al.ll2005l ). and from A 0620-00 (|Gallo et all 120061 ) 
providing support to the idea that an outflow also plays a 
significant role in quiescence. 

To facilitate comparison between developing models 
and data we present here the most complete SED available 
for V404 Cyg. It is primarily based on a simultaneous cam- 
paign performed in 2003 using Chandra, HST, and several 
optical and radio observatories. We supplement the cover- 
age with non-simult aneous optical and near-IR data from 
Casares et al.l Il993l) and near to mid-IR Spitzer data from 
Muno fc Mauerhanl l|2006h and lGallo et~aD (120071). Previous 
result s from th is campaign were p resented by Hynes et alj 
l|2004h : see also lCorbel et al] l|200Sl ). 

2 OBSERVATIONS 
2.1 HST ultraviolet data 

Ultraviolet photometry used the Advanced Camera for Sur- 
veys fACS: lBofn et alfcOOTl ) onboard the Hubble Space Tele- 
scope (HST). The observations performed are summarised in 
Table fT] V404 Cyg is a faint UV source and reddened, so five 
orbits were used covering the F250W and F330W filters. 

We used standard ACS techniques to reduce the data. 
Individual images had been pre-processed to the flat-fielded 



stage using the automatic pipeline, CALACS; we found no 
need to refine this reduction. We then combined all of the 
images offline into a geometrically correcte d master image 
using Multidrizzle |Koekemoer et all [2002) with standard 
settings. We found excellent registration between individual 
images, so no need to improve the shifts applied. 

We performed photometry using the IDL/AstroLib 
aperture photometry routine APER. As V404 Cyg is faint 
we used a relatively small aperture to perform photometry, 
of radius 0.125", with background defined by a 2.5" annu- 
lus. No stars in the field were bright enough to determine 
per-field aperture corrections , so we use the tabulated val- 
ues from lSirianni et al] (|2005l ). This aperture collects about 
75 % of the total light, and almost all of the sharp core of 
the point spread function. We measured the position of the 
target in the F330W bandpass and then fixed it for F250W, 
after verifying that it was consistently placed in the two. 
We note in passing that we do detect both V404 Cyg and 
the nearby contaminating star. We find an offset of the con- 
taminating star relative to V404 Cyg of Aa = +0.005 sec 
(= 0.06") and AS = 1.43", an improvement on the rela- 
tive astrometry of t he stars compared to that reported by 
ICasares et ail (|l993l ). 

A significant source of systematic uncertainty in 
ACS photometry is charge transfer ineffici ency (CTI) due 
to radiation damage to the CCDs (jRiess fc Mackl 120041 : 
iPavlovskv et all 120051) . although fortunately the observa- 
tions were performed early enough in the lifetime of ACS 
that CTI is unlikely to be a severe problem. The effect is 
worst at low light levels, so maximised for faint UV sources. 
CTI is also least well calibrated f or these cases, and the 
correction prescriptions provided bv IPavlovskv et al.l (|2005l ) 
formally diverge for negligible sky background or source 
counts. We face both problems with these data. Sky back- 
grounds in the near-UV are extremely low, dominated by 
zodiacal light and earthshine in the F330W bandpass, and 
geo-coronal emission at shorter wavelengths. Both were es- 
sentially unmeasurable in our images, and the background 
is expected to be dominated by CCD dark current, pro- 
ducing 0.5-2.0 e~ /pixel for our exposures. We attempted 
to estimate the CTI losse s in two ways. Firstly , we fol- 
lowed the prescription in IPavlovskv et al] l|2005l ). except 
that we used the estimated background (including the dom- 
inant dark current term) from the ACS Exposure Time 
Calculator. This prescription predic ted losses of ~ 5 % . 
Based on the actual data shown by iRiess fc Mackl (|2004l ) 
however, this prescription appears to overestimate the CTI 
losses for the faintest sources. For the low background case 
(0.5-1.5 e~/pi xel) at source brightn esses up to a few hun- 
dred electrons, IRiess fc~M ack ( 2004]) found losses of around 
0.035 mag at maximum distance from the readout amplifier. 
If we rescale this to the times of our observations and po- 
sitions of the source we expect losses of about 2.5 %. We 
use these latter estimates to correct our measured source 
brightnesses and assign an additional error estimate of 2.5 % 
to account for the uncertainty in applying CTI corrections. 
The uncertainties in the CTI corrections are much smaller 
than the statistical uncertainties in the measured fluxes, so 
this is not the dominant source of uncertainty. 

We tabulate our photometric results in Table [5] Since 
we are measuring counts across a broad bandpass rather 
than monochromatic fluxes, we tabulate the data in two 
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Table 2. Results of HST photometry 



Band Count rate 



Flux 

(erg cm" 2 s" 1 A" 1 ) 



Observed 

F250W 0.109 ±0.042 (0.52 ± 0.20) x 10~ 18 
F330W 1.06 ±0.14 (2.37 ± 0.30) x 10" 18 

Synthetic KOm 

F250W 0.073 

F330W 0.57 

Synthetic KOm, truncated at 4200 A 

F250W 0.024 

F330W 0.57 



ways; the corrected number of electrons per second, and the 
average flux per unit wavelength using photometric zero- 
points estimated by the ACS pipeline. All have been cor- 
rected for CTI, and to a nominal infinite aperture. 



2.2 IAC80 photometry 

We obtained simultaneous R band observations using the 
Thomson CCD camera at the 0.8 m IAC80 from 20:50-05:51 
UT on the night of 2003 July 28/29. We used 300 s exposures 
with 52 s readout time between them and bias-corrected and 
flat-fielded all images in the standard way using iraeE 
Seeing conditions were not good enough to cleanly sepa- 
rate the contributi on of the target and it s nearby (1.4") 
line-of-sight star l|Udalski fc KaluznH Il99l[ ) so we applied 
straightforward aperture photometry using a large aperture 
of 3.5"which adds the flux from both stars as described 
by IZurita et al] \2004) . We determine an average magni- 
tude R = 16. 6 3 ± .01 in acceptably close agreement with 
ICasares et ail (|l993T ). indicating little difference in the opti- 
cal brightness of the system compared to earlier data. 



2.3 WHT optical spectroscopy 

To reconstruct the optical SED of the source we use the 
WHT data from the program. This covered a larger range 
of wavelengths than the Gemini data, and was simultane- 
ous with the HST observations. These observations used the 
ISIS dual-arm spectrograph on the WHT. To maximise effi- 
ciency and minimise readout time and noise, we used the sin- 
gle red-arm mode with the R316R grating and MARCONI2 
CCD. We set exposure times to 200 s, with ~ 17 s dead-time 
between exposures. To maximise photometric accuracy, we 
used a 4" slit, so our spectral resolution was determined 
by the seeing (median ~ 1.3"), and was typically ~ 5.5 A 
(250 km s" 1 ). We performed bias correction and flat field- 
ing using standard IRAF techniques. The slit was oriented 
in the same was as in our previous observations and cov- 
ered the same comparison star. We extracted spectra of both 
of these stars, and the nearby blended star, with the same 

1 IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



techn iques previously described (|Hvnesll2002l ; IHv nes et al.l 
2002). We performed wavelength calibration relative to a 
single observation of a CuNe/CuAr lamp and corrected for 
time-dependent variations in the wavelength calibration us- 
ing Telluric absorption features. We calibrated the on-slit 
comparison star relative to Kopff 27 (|Stonelll977h . and then 
calibrated all spectra of V404 Cyg relative to this on-slit 
comparison. 

To compare the WHT data with that from the IAC80 
we perform synthetic photometry convolving the spectrum 
with a R bandpass. We found an offset between the two, 
most likely due to systematic errors in the WHT calibration. 
Since the IAC80 photometry agrees well with earlier mea- 
sures we assume that this calibration is correct and rescaled 
the WHT data up by a factor of 1.3 to agree with the IAC80 
R band data. 



2.4 Gemini data 

We used additional optical spectroscopy obtained with the 
R831 grating and standard EEV CCDs of the CMOS spec- 
trograph on Gemini-N to extend our continuum lightcurve. 
Exposure times were 40 s and we binned and windowed the 
images to reduce the dead-time between exposures to 12 s. 
With a 5" slit and 1.1" median seeing we obtained a spectral 
resolution of 5.0 A (230 kms -1 ). We performed data reduc- 
tion, spectral extraction, and wavelength and flux calibra- 
tion in the same way as for the WHT data. Wavelength 
calibration used a CuAr lamp, and we performed flux cali- 
bration relative to the same on-slit comparison star as used 
for the WHT observations, and applied the same rescaling 
to ensure consistency with IAC80 photometry. 



2.5 Chandra X-ray data 

We use data taken from two recent Chandra o bservations of 
V404 Cyg, one on 2000 April 26 for 10,295 s l|Garcia et all 
l200ll . iKong et al.ll2002l). and another o n 2003 July 28/29 
for a total of 61.200 s (|Hvnes et al.ll2004l ). The first observa- 
tion used 1/4 sub-array mode, with time resolution of 1.14 s; 
while the second used 1/8 sub-array mode, with a resolution 
of 0.4 s. In the original set, 1587 counts were reported, while 
the second had 1941 counts. We reanalyzed the standard 
pipeline-processed level-2 data from both data sets in CIAO 
v3.2. The extraction aperture was set to a radius of 6-pixels 
and we retained only events between 0.3-7.0 keV to reduce 
the background. The background can be particularly signif- 
icant for the S3 chip, so we used large background regions 
with a 48-pixel radius. The background produced approxi- 
mately 4.6 counts in the source aperture for the 2003 obser- 
vation, and 0.2 counts in the 2000 observation. 



2.6 VLA radio data 

VLA data were obtained in the A configuration from 29 July 
00:20 UT until 14:15 UT. The data were reduced and each 
time bin imaged using standard procedures within AIPS. 
Absolute flux calibration was obtained relative to 3C 286. 
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2.7 WSRT radio data 

The Westerbork Synthesis Radio Telescope (WSRT) ob- 
served V404 Cyg from 28 July 17:13 UT till 29 July 
05:33 UT, corresp onding to the period D1-F1-D2-F2 in 
iHvnes et al] {2004) and labelled in Fig. [1] We used two fre- 
quencies: 4.8 and 1.4 GHz, with ~ 5 hours on source per 
frequency. The first hour was unusable because the source 
was too low on the horizon. The full 4.8 GHz observation 
(with useful coverage 18:40-05:13 UT) gives a flux density 
of Stot = 0.19 ± 0.03 mjy. At 1.4 GHz the average flux den- 
sity was Stot, 1.4 = 0.18 ± 0.08 mjy. 



2.8 Published Non-simultaneous Data 

We supplement our dataset with non- simultaneous points 
from the literature to fill in gaps in coverage. The most 
comprehensive cove rage of the optical-IR SED is from 
ICasares etahl (|l993h . who present BVR and JHK photom- 
etry from a coordinated campaign obtained on 1990 June 10. 
These points do raise some concern that they were taken 
soon after the outburst and may not be fully quiescent. 
There is also a large gap, 13 years, relative to the recent 
data. Nonetheless, we find excellent agreement in R between 
our simultaneous IAC80 photometry and the earlier data. 
Re cent Spitzer data in the m id-IR have also be en presented 
bvlMuno fc Mauerhanl {2006 ) and reanalysed bv lGallo et al.l 
(2007). We include these 4.5, 8.0, and 24 /zm photometry in 
the SED using Gallo's data which differ significantly only in 
the 24 /im measurement. 



IAC80 R band 



3 THE UV, OPTICAL, AND IR SPECTRAL 
ENERGY DISTRIBUTION 

3.1 Contamination of the R band by Ha 

Using the WHT spectrum we can estimate the contribution 
of line emission to the photometry and correct for this effect. 
The strongest line is Ha with an equivalent width of 16 A. 
This translates to a change in R magnitude of 0.01 mag, and 
we correct for this befor e using these data in the SED. 

ICasares et al.l (|l993l ) found a stronger line during 1990- 
91 with average equivalent width 38.7 A. This suggests a 
larger contribution to the earlier photometry of ~ 2 %, and 
we use this value to correct the photometry of lCasares et al.l 
(1993). We note that in both cases this contamination is 
effectively negligible in comparison to other uncertainties. 



3.2 Constraints on the optical disc contribution 
from variability 

We expect the optical emission in quiescent BHXRTs to be 
a combination of light from the disc and light from the com- 
panion star. In the case of V404 Cyg, this is overwhelm- 
ingly dominated by the companion. Traditionally measures 
of the disc fraction use veiling of the photospheric lines of 
the companion star and are dependent on having high res- 
olution spectra at the time of the observation. They can be 
rather sensitive to assumptions abou t the spectrum of the 
companion star (Hynes et al. 2006). ICasares et al.l (|l993l ) 
investigated several methods of fitting the veiling at differ- 
ent wavelengths within the optical spectrum. Despite large 
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Figure 1. Optical, X-ray, and radio lightcurves of V404 Cyg dur- 
ing our simultaneous observations. The upper panels show the full 
flux observed by IAC80 photometry and WHT and Gemini con- 
tinuum spectroscopy. The dashed line in the second panel is a lin- 
ear fit to the lower-envelope. The third panel shows the WHT and 
Gemini continuum data after subtraction of this lower-envelope. 
The photometric data has been treated in the same way after 
re- normalising to agree with the WHT data and is shown dotted. 
The fourth panel shows the X-ray data and the fifth the radio 
lightcurve. 



uncertainties in their results, they conclude on approximate 
values of ra = 0.26, rv = 0.15, and tr = 0.09, where the r 
values refer to the frac tional contribut i on of the disc to the 
totafl In the near-IR. IShahbaz et al] (1996) put an upper 
limit on the K band veiling of tk < 0.14. 

We cannot make a reliable determination of the disc 
fraction from our low resolution spectroscopy, but can in- 
stead use a different method to place a lower limit on the 
disc fraction. The optical continuum is significantly vari- 
able, and the variability presumably must be attributed to 
the disc. If we assume that the optical continuum lightcurve 
consists of three components: companion star, slowly vary- 




2 ICasares etHl j l993h quote slightly different values as they de- 
fine the veiling as the disc light as a fraction of the stellar light 
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ing disc light, and nickering disc light, then the fraction of 
light due to the flickering component provides a robust lower 
limit on the disc fraction. 

We show IAC80 photometry together with WHT and 
Gemini continuum lightcurves in Fig. [T] Our observations 
were centre d around photometric p hase 0.65 using the 
ephemeris of lCasares fc Charles! (|l994h implying that the el- 
lipsoidal modulation should be rising from the deeper mini- 
mum at phase 0.5 to the maximum at phase 0.75. This rising 
trend and the approximate amplitude of ellipsoidal modula- 
tions expected are consistent with the general upward trend 
we observe. For a long period system such as V404 Cyg, we 
can approximate the ellipsoidal modulation within a night 
by a linear trend so we remove the slowly varying light by 
subtracting a linear fit to the lower envelope of the contin- 
uum lightcurve. We find that the remaining flickering com- 
ponent is about 4 percent of the total light in both the WHT 
and Gemini segments of the lightcurve and correlates very 
well with t he X-ray variation s (as also found for the Ha 
emission by iHvnes et al.|[2004T ). This 4% flickering light is 
actually not far off of spectroscopic estimates based on veil- 
ing near Ha as discussed above. This makes it plausible that 
flickering light source does contribute most of the disc light 
in the optical. 




14.2 14.4 14.6 14.8 

Log v (Hz) 

Figur e 2. Optical to near-IR SED using photometry from 
ICasares et al.l lll993l) . No dereddening has been applied to the 
data. The solid line is the best fit of a K0 HI spectrum reddened 
with Ay = 4.04. The dotted line corresponds to Ay = 3.6 and 
the dashed line to Ay = 4.4. The solid lines below the BVR 
and K points are the ranges of plaus ible companion con tribu- 
tions to the SED based on estimates o f C asa.res et al.l (Il993h and 
IShahbaz et ah! dl996r i. 



3.3 The stellar contribution 

We do not have sufficient information to independently 
constrain the stellar contribution to the SED at each 
wavelength, so opt to instead model the stellar spec- 
trum based on est i mates of the system parameters. From 
ICasares fc Charles! (I1994T) we have a spectral type estimate 
of G8-K2III-V. IShahbaz et all (jl994l ) estimates companion 
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star parameters of A/2 = O.Ttg'a M© and R2 
implying logg = 2.76 ± 0.08 and a giant classification. The 
SED of a giant star should be the best approximati on for 
V404 Cyg, although strictly it is a stripped giant (|Kind 
1 19931 ). For a KOni spectral type we find T eB ~ 4570 K u sing 
the temperature scale derived by Ivan Belle et al.l (|l999T l for 
G, K, and M giants. We model the spectrum of the com- 
panion by interpolating on t he grid of model atmos phere 
spectra of late-type giants of lHauschildt et al.l l|l999l ) after 
loga rithmically re-binnin g the models to Alogv = 0.01. 

ICasares et al.l l|l993h estimated the interstellar extinc- 
tion to be Ay = 4.0 based primarily on comparing the 
(B — V) colour with that expected for a K0 ill star after cor- 
recting for their estimates of disc contamination. We tested 
this by fitting the model K0 ill spectrum described above 
to the full set of BVRJHK magnitudes (folding the model 
spectrum with bandpasses to generate synthetic photome- 
try) constraining the fit to preferably lie wit hin ranges im- 
plied b y th e disc contamination constraints of lCasares et al.l 
\ 19931 ) and lShahbaz" et al.l (1 19961 ) and shown in Fig. [5] by ver- 
tical bars. Our fit confirms the earlier estimate with a formal 
best fit for Ay = 4.04, shown in Fig. [2] The residuals scatter 
both above and below the line indicating that intrinsic vari- 
ability and/or calibration uncertainties are playing a role, so 
we do not attempt to estimate an error on the extinction but 
we do show companion spectra reddened by Ay = 3.6 and 
Ay — 4.4 to indicate the sensitivity to extinction. Extinc- 
tions much outside of this range would require attributing 



a larger fraction of the light to the disc than has previously 
been estimated. 

The alternate G8 and K2 possibilities are almost indis- 
tinguishable from K0, so we do not show them. The KOni 
stellar spectrum clearly provides a good description of the 
UV to near-IR data, consistent with previous estimates that 
the companion dominates the light in both the optical and 
infrared. We expand the view in Fig. |3]to include both near 
to mid-IR d ata from Spitzer and near- UV data from HST. 
As found bv lMuno fc Mauerha n (2006), an IR excess does 
appear to be seen by Spitzer. A weak UV component also 
seems to be present, although thi s is not as pronounced as 
seen i n other quiescent BHXR Ts (|McClintock fc Remillardl 
120001 : iMcClintock et al.1 120031 : Hynes et al. in prep.) The 
F250W data (the highest frequency) shown in Fig. [3] has 
already been corrected for an estimate of the red-leak as 
will be discussed in Section f3. 5 1 



3.4 The mid-IR excess 

The Spitzer data clearly lie above a plausible extrapolation 
of the companion star spectrum alone (Fig. [3}. Given the 
insensitivity of the IR spectrum to the temperature of the 
companion and interstellar extinction, we cannot adjust the 
companion spectrum in any way to fit these IR data. While 
these were not simultaneous with the optical data, we do 
not expect the continuum brightness of the companion to 
vary from epoch to epoch. Small variations due to ellipsoidal 
modulations will be present as a function of orbita l phase, 
but these are less than ±10 % JShahbaz et al. I [19941 ). 

uno fc Mauerhanl (|2006h interpreted the mid-IR ex- 
cess seen in V404 Cyg as the tail of a Shakura-Sunyaev 
accretion disc. This model, however, predicts a UV excess 
which is inconsistent with our HST observations, exceeding 
our F330W measurement by an order of magnitude and the 
F250W by two orders. During quiescence instead most of the 
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Figure 4. Power-law fits to Chandra data from 2000 (left) and 2003 (right). 




13.0 13.5 14.0 14.5 15.0 

Log v (Hz) 

Figure 3. UV-Optical-IR SED. No dereddening has been applied 
to the data. The short solid line segment and filled symbols in- 
dicate new data, open symbols mark archival photometry. The 
highest frequency HST point has been approximately corrected 
for red- leak (Section [33}. The dashed line is a model K0III SED 
reddened with Ay = 4.04. The dotted line is a 2000 K blackbody 
component adjusted so that the sum of the two (shown with a 
long solid line) fits the short wavelength Spitzer data. 

disc is expected to be in a cool state, 2000-30 00 K, with a 
more isothermal temperature distribution (see iHvnes et al.l 
120051 and references therein). The contamination from this 
material will then be most pronounced in the mid-IR as it 
is cooler than the companion star, but no UV excess will be 
present. 

We use a solid line in Fig. [3] to show a model includ- 
ing a 2000 K blackbody in addition to the companion star 
spectrum discussed above. This can adequately fit both the 
near-IR to UV SED and the Spitzer near-IR data (it is less 
successful at 24 fim). The blackbody flux is consistent with 
a di sc of radius ~ 9 x 1 11 cm tilted at an inclination of 
56° ijShahbaz et al.l ll994) and at a distance of 3.5 kpc. For 
the system parameters discussed above, this corresponds to 
about 60 percent of the radius of the Roche lobe of the black 
hole. This is a very plausible size for a quiescent accretion 
disc. 



We note that while lMuno fc Mauerhanl (|2006l ) preferred 
a Shakura-Sunyaev disc model for the IR excess in V404 Cyg 
specifically, in other systems they invoked a much cooler 
circum-binary disc. Such a disc would remain a possibility 
in the light of our UV observations as it would not contribute 
to the UV flux. 

The difficulty in fitting the 24 /im point supports a 
more controversial interpretation as proposed bv lGallo et al.l 
l|2007l ). The flux in this band is rather close to an extrapo- 
lation of the flat radio spectrum raising the possibility that 
the IR excess could be due to synchrotron jet emission. The 
major difficulty in quantifying this is that the radio is known 
to be strongly variable (Section [SJ and was not observed si- 
multaneously with the 24 /^m point. We will return to this 
topic in Section \6. 31 



3.5 The UV excess 

Our data also provide tenuous evidence for a UV excess. 
One possibl e explanation for this is r ed-leaks in the UV 
bandpasses (|Chiaberge fc Siri anni 2007). Although this may 
make a large contribution in F250W, the effect is not ex- 
pected for the F330W filter. The estimates for red-leaks from 
IChiaberge fc Siriannil l|2007l ) are of no use as they only con- 
sider unreddened stars, so we perform our own calculations. 
In addition to red-leaks, because the spectrum is so steep in 
the UV it is hard to gauge the expected flux in a broad-band 
without performing synthetic photometry. 

To better assess the significance of any UV excess, 
we therefore use synphot (jLaidler et al.l 120051 ) to estimate 
count-rates in each bandpass for our model reddened KOm 
spectrum. We consider two spectra, one with the full spectral 
coverage, and the other truncated above 4200 A to disentan- 
gle in-band stellar flux from the red-leak. We find that the 
red-leak is negligible in the F330W bandpass, as expected, 
and that the in-band stellar flux can account for approxi- 
mately half of the count-rate seen in this band; this is a little 
more than a 3 a excess. The F250W band appears to have 
a large contribution from red-leak, about half of the flux in 
the band. After subtraction of the companion star, both in- 
band and red-leak, the residual excess is 0.36 ± 0.40 s _1 , so 
is not a significant detection of disc light in this band. 
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We therefore appear to have a significant detection of a 
UV excess in the F330W band only, but even this is at the 
mercy of uncertainties in the companion spectral classifica- 
tion (which do become a concern in the UV) and reddening. 
We therefore choose not to pursue the spectrum of the UV 
excess further. 



4 THE X-RAY SPECTRA AND VARIABILITY 
4.1 X-ray Spectra Analysis 

We extracted spectra for both Chandra observations with 
CIAO v3.2 and analysed them with XSPEC vl2.2. We used 
standard CIAO techniques to select response files accord- 
ing to the CCD temperature. We considered many spectral 
models within XSPEC, all with interstellar absorption, and 
report results from fitting with power-law, bremsstrahlung, 
and Raymond-Smith models here. Blackbody fits to the data 
were so poor that we neglect them for discussion. Due to the 
lower count rate in the 2003 observation, we were prompted 
to try using CASH statistics in tandem with the Churazov 
approximation for chi-squared statistics (Cash, 1979; Chura- 
zov, 1996). For each set of statistics to be used, we grouped 
the spectra into at least 15 counts per spectral bin; further 
binning results in a loss of spectral information. Churazov 
weighting estimates the weight for a given channel by aver- 
aging the counts in surrounding channels. CASH statistics is 
a method designed to estimate the best-fit parameters using 
unbinned or slightly binned data, which can be particularly 
useful when the source yields few photons. We also checked 
the consistency of CASH with unbinned data, and found the 
results to be equivalent. 

All models except the blackbody model give acceptable 
fits to both data sets with the power-law providing the best 
fit. Both statistical methods provide res ults that are equiv- 
alent to previously published results bv lKong et alj (|2002l ) 
for the 2000 observation. The best-fitting model is shown in 
Fig. HI For t h e 200 3 observation, we find similar results to 
ICorbel et al.l |2008), and a plot is also displayed in Fig. [4] 
Table |3] summarises the results of the spectral fitting. 



4.2 Color-Intensity Diagram 

To test for any variation in the spectrum with luminos- 
ity we construct a hardness-intensity diagram. We sepa- 
rated the data into two bandpasses of approximately equal 
count rates, which were 0.3-1.75 keV (soft) and 1.75-7.0 keV 
(hard) for each set. To account for the reduction in sensi- 
tivity between the two epochs, we used an assumed fixed 
model of a photo-absorbed power-law to calculate the ra- 
tio of expected count rates in these bandpasses between the 
two observations and hence rescale the 2000 observations to 
match those from 2003. The resulting hardness-intensity di- 
agram is shown in Fig. [5] We see no variation in hardness 
either within or between the two observations, the first of 
which had a substantially higher count rate. This indicates 
no detectable spectral change over a factor of ten in source 
luminosity. 
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Figure 5. Hardness— intensity diagram from both Chandra ob- 
servations. The two datasets were combined using the power- 
law model to estimate the change in sensitivity between the two 
epochs. Solid circles indicate 2000 data and open circles that from 
2003. 



5 RADIO SPECTRA AND VARIABILITY 

To constrain the radio spectral shape we first examine each 
dataset individually. We obtain a = —0.07 ± 0.16 (where 
/„ oc v a ) from the VLA data and a = +0.04 ± 0.38 from 
the WSRT. The difference between them is not statistically 
significant so we also perform a joint power-law fit to both 
datasets, allowing the relative normalisation of the two to 
vary as they sample different portions of the lightcurve. Pri- 
marily this is constrained by the VLA data as the uncer- 
tainties in the WSRT fluxes are larger. We find a spectral 
index of a = —0.05 ± 0.15. All the fits are consistent with 
the canonical flat spectrum as commonly seen in hard state 
black hole binaries. 

The VLA data exhibit substantial variability. We show 
the higher quality 8.5 GHz lightcurve together with that 
from Chandra in Fig.[T] The radio is clearly strongly variable 
and shows both large flares and dips not dissimilar to those 
seen by Chandra. The most dramatic feature rises from al- 
most zero flux to the highest peak in ab out 30min, as was 
also th e case in the dataset presented by iMiller- Jones et al.l 
(2008). No clear correlation between radio and X-ray is ap- 
parent to the eye, however. To further test this, we calculate 
the cross-correlation function between the X-ray and radio 
data and show this in Fig. [S] No compelling correlation (or 
anti-correlation) is present although an unremarkable peak 
is present at zero lag. 

Since WSRT is a linear array, it covers a full synthesis in 
12 hours and hence is not capable of the high time-resolution 
achieved by the VLA. Nevertheless the source does show sig- 
nificant variability over this period, even if a possible cor- 
relation with X-ray /optical is difficult to infer. When the 
observation is split in 2 sub-intervals, 18:40-24:00 and 00:00- 
05:13, the resulting flux densities are: Si = 0.25 ± 0.04 mjy 
and £2 = 0.14 ± 0.04 mjy respectively. The second interval, 
covering D2-F2, overlaps with the first part of the VLA ob- 
servation, which shows a significant increase at 8.5 GHz be- 
ginning at 03:00 UT (with a peak flux density of 0.6 mjy). 
However, the VLA flare only lasts for ~ 1 hour, of which only 
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Table 3. Best-fitting parameters for Churazov fits to Chandra data. 



Model N H a kT x 2 /d-o.f. 0.3-7.0 keV Flux 

(10 22 cm- 2 ) (kcV) (ergem^s- 1 ) 
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Figure 6. Cross correlation function of radio lightcurve with re- 
spect to the X-rays. Positive lags would indicate radio lagging 
behind the X-rays. Dotted and dashed lines are 1 — a and 3 — u 
expectations for uncorrelated variations. 



about one-half was covered by WSRT. Finer sub-intervals 
proved impractical as there was insufficient coverage of the 
UV plane to isolate flux from V404 Cyg. 



6 THE BROAD BAND SPECTRAL ENERGY 
DISTRIBUTION 

6.1 Comparison with existing data 

Previous sections of this paper have focused primarily on 
one or two bands; we will broaden our scope here to ques- 
tions best addressed with the whole SED. We show the com- 
plete radio-X-ray SED in Fig. [7] We have also compiled 
all the data points, as plotted, in Table [4] to allow easier 
comparison with future models. We include a point repre- 
senting just the nickering component in the optical to rep- 
resent a lower limit on the acc r etion light. We adopt the 
exti nction curve of Fitzpatrick (1999) in the optical and 
UV. Ilndebetouw et al. | (|2005h for the Syitzer/lKkC d at a and 
assume Ai^mj-^K — 0.5 following IChapman et alj f2 009b 
Our d ataset should supersede that used by iNaravan et al.l 
4 1997ft for future modelling of V404 Cyg. Compared to ear- 
lier work we now have simultaneous X-ray and optical data, 



shorter wavelength coverage extending into the near-UV, 
and simultaneous radio observations. 

In Fig. [5] we illustrate th e additional constrain ts now 
possible, showing Model 1 of INaravan et al.l (|l997| ) over- 
laid on our SED after scaling downwards to fit our X-ray 
spectrum. This is of course not intended to be a detailed 
model fit. Not only should a new model be computed to fit 
a lower X-ray luminosity, but there have been refinements in 
modelling of advective accretion flows in the last decade as 
well. Our main point is to stress the strong constraint pro- 
vided by the low near-UV flux. Model 1 was an accep table 
fit to the SED data available to INaravan et all (|l997ft but 
clearly over-predicts the near-UV substantially compared to 
our new data. Alternative m odels with winds presented by 
iQuataert fc Naravanl (1 19991 ). on the other hand, would be 
consistent with the UV data. This constraint on the UV to 
X-ray flux ratio can be used to constrain current and future 
models. 



6.2 The origin of optical flickering 

The flux of the optical flickering component falls at an in- 
triguing location in the SED (Fig. [TJl . It is a little above a 
straight extrapolation of either the X-ray or radio power- 
laws. An extrapolation of the radio is plausible, as a flat- 
spectrum is expected to continue until a break to optically 
thin synchrotron at higher frequencies. We would not ex- 
pect a straight extrapolation of the X-ray s pectrum at all, 
and a ll models for the broad-band SED (e.g. INaravan et al.l 
Il997ft involve curvature between X-ray and optical bands.) 
Here correlations are more telling evidence than fluxes. We 
find an excellent correlation between X-ray and optical con- 
tinuum (Fig. [TJ, but no clear correlation between X-ray 
and radio (and by extension optical). This suggests that 
the optical variability and the X-ray emission are closely 
connected, either by direct emission from the same region, 
as expected from early advection dominated accretion flow 
(ADAF) models, or indirectly through X-ray irradiation o f 
the disc, as inferred for Ha emission bv lHvnes et al.l l|2004ft . 
The X-ray and radio emission, on the other hand, do not 
obviously appear correlated on the timescales observed (al- 
though a correlation on longer timescales is certainly pos- 
sible) suggesting that the mechanisms for variability in the 
two are not connected. 
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Figure 7. Broad band SED. Solid points are simultaneous data from our campaign, open points are non-simultaneous. Triangles are 
used for radio data, squares for Spitzer, circles for optical photometry, and stars for HST UV photometry. Solid lines indicate the stellar 
fit to the optical data and the power-law fit to the radio. 



6.3 The extent of synchrotron emission 



There is clearly a mid-IR excess prominent in the 24 fim 
band. While it is possible to partially explain this with emis- 
sion from a cool disc, synchrotron emission fro m the radio 
jet p rovides a better description of the excess dGallo et al.l 
120071 ). An extrapolation of the flat spectrum beyond the IR 
region is problematic for several reasons. We have argued 
above that the optical variability appears directly or indi- 
rectly associated with the X-ray emitting component rather 
than with the radio. In addition, the low UV flux observed 
is inconsistent with a straight extrapolation of the flat spec- 
trum all the way to the UV, and problematic even if the 
spectrum is assumed to break in the optical. This argues 
that a break in the jet spectrum at longer wavelengths is 
needed so that the non-stellar optical/UV emission origi- 
nates from the disc rather than the jet. The most likely 
interpretation is then that the break from fiat spectrum to 
optically thin synchrotron occurs in the mid to far-IR, al- 
though given the large amplitude of radio variability only 
simultaneous radio and mid-IR observations could address 
this question conclusively. 



7 CONCLUSIONS 

We have compiled the comprehensive dataset describing the 
spectral energy distribution of the black hole X-ray binary 
V404 Cyg in quiescence. This includes new and simultane- 
ous radio, optical, UV, and X-ray data, supplemented with 
additional non- simultaneous data in the optical through in- 
frared. 

Our specific findings are: 

(i) There is no need for a model to include a substantial 
contribution from accretion light over the near-UV to opti- 
cal spectral range. The SED over this range (including the 
vacuum UV) is well described by a K0 IV star reddened by 
Ay = 4.0. The very low UV flux detected only marginally 
supports an accretion component, and rather tightly con- 
strains some spectral models of quiescent accretion flows. 

(ii) As suggested bv lMuno fc Mauerhanl (|2006l ) no plausi- 
ble companion spectrum and interstellar extinction can fully 
explai n the mid-IR; an IR excess from a cool disc or possibly 
a jet (|Gallo et alj|2007h appears to be required. 

(iii) The X-ray spectrum is consistent with a T ~ 2 power- 
law as found by all other studies to date. There is no evi- 
dence for any spectral variation over the range of a factor of 



10 R. I. Hynes et al. 



-10p 




-14 r 



-15 1 L ■ i i i i i i 

12 13 14 15 16 17 18 

Log V (Hz) 

Figur e 8. Comparison of our SED with Model 1 of lNaravan et all 
1119971 ). The dashed line is their ADAF model, the dotted line is 
now just the stellar emission. The solid line is the composite of 
the two. 

10 in luminosity that has been observed over two epochs by 
Chandra. 

(iv) The radio shows a flat spectrum in f v as is commonly 
seen in hard state X-ray binaries and associated with a com- 
pact jet. The break frequency between a flat and optically 
thin spectrum most likely occurs in the sub-mm to far-IR 
regime, but is not conclusively constrained by existing data. 
The radio is substantially variable but puzzlingly there is no 
clear correlation with X-ray variability. If the radio variabil- 
ity reflects changes in the jet power, this suggests that the 
X-ray variations occur independently of changes in the jet, 
or with a very large lag (many hours) between the two. 

Overall, the likely interpretation of the SED is as fol- 
lows: 

(i) The KOiii companion star dominates the light from 
near-IR to near-UV wavelengths. 

(ii) The cool outer regions of the accretion disc 
may produce a modest mid-IR excess seen by Spitzer 
i|Muno fc Mauerhan|[200l) . 

(iii) The inner accretion flow produces highly variable X- 
ray emission but apparently only weak UV emission. This 
irradiates the a ccretion disc and is responsible for emission 
line variability l|Hvnes et alj 120041, The optical continuum 
variability may also arise from reprocessing but direct emis- 
sion from the inner accretion flow remains a possibility. 

(iv) A compact jet is driven from the inner accretion disc 
producing variable flat-spectrum radio emission likely ex- 
tending to the far and possibly mid-IR. The size of the jet is 
smaller than 4.5 AU (about 30 x the binary s eparation) for 
a distance of 3.5 kpc ({Miller- Jones et alj|200a ). 
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Table 4. Compilation of SED data used in Fig. [7] 
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